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Abstract 

In this lecture, the early history of asymptotic freedom is discussed. The first 
completely correct derivation of (3$ in non-abelian gauge theory (Khriplovich, 1969) 
was done in the Coulomb gauge; this derivation is reproduced (in modernized terms) 
in Sect. [2 A qualitative physical explanation of asymptotic freedom via chromomag- 
netic properties of vacuum (Nielsen, 1981) is discussed in Sect. El 

1 Introduction 

Evolution of the QED coupling a(/x) is determined by the renormalization group equation 

= ~ 2 ^ a ^ a ^ ' ft** = ^ + & (£) 2 + • ■ ■ ^ 

where the 1-loop /3-function coefficient is 

0o = ~ (1-2) 

(we use the MS renormalization scheme). The sign (3q < corresponds to the natural 
picture of charge screening — when the distance grows (/x reduces), a(/x) becomes smaller. 

Non-abelian gauge theories have (3q > 0. This corresponds to the opposite behaviour 
called asymptotic freedom — the coupling becomes small at large /x (small distances). 
Discovery of this fact came as a great surprise, and completely changed ideas about ap- 
plicability of quantum field theory for describing the nature. In fact, it was done several 
times, by several researchers, independently. 

The first paper in which (3$ > was obtained is by Vanyashin, Terentev [T]. They 
considered non-abelian gauge theory with the gauge boson mass term introduced by hand. 
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Of course, this theory is non-renormalizable. It contains longitudinal gauge bosons whose 
interaction grows with energy. They obtained 



The contribution — 1/6 comes from longitudinal gauge bosons. In the massless gauge 
theory it is canceled by the ghost loop (which was not known in 1965). Now we know that 
renormalizable massive gauge theory can be constructed using the Higgs mechanism. In 
this theory, the result (II .3p is correct; the contribution —1/6 comes from the Higgs loop. 

The first completely correct derivation of /3q in non-abelian gauge theory has been 
published by Khriplovich [2j. He used the Coulomb gauge which is ghost-free. The Ward 
identities in this gauge are simple (as in QED), and it is sufficient to renormalize the 
gluon propagator in order to renormalize the charge. This derivation clearly shows why 
screening is the natural behaviour, and how non-abelian gauge theories manage to violate 
these general arguments. It is also very simple. We shall discuss this derivation of Po in 
Sect. El following [2J, but using a modernized language. 

't Hooft discovered asymptotic freedom in 1971-72 while studying renormalization of 
various field theories in dimensional regularization [3]. He reported his result during a 
question session after Symanzik's talk at a small conference in Marseilles in June 1972. 

And finally, Gross, Wilczek [I] and Politzer [5] discovered it again in 1973. They were 
first who suggested to apply asymptotically free gauge theory (QCD) to strong interactions, 
in particular, to deep inelastic scattering. 

The early history of asymptotic freedom is discussed in several papers, e.g., O [7J [8], 
and in the Nobel lectures [9]. 

Now Pq is derived in every quantum field theory textbook. In the standard approach, 
the covariant gauge is used (with Faddeev-Popov ghosts, of course). The coupling constant 
renormalization can be obtained from renormalizing any vertex in the theory together with 
all propagators attached to this vertex. Usually, the light-quark - gluon vertex or the ghost 
- gluon vertex is used (calculations are slightly shorter in the later case). The 3-gluon vertex 
or even the 4-gluon one [10] can also be used. All these standard derivations can be found, 
e.g., in [UJ. The infinitely-heavy-quark - gluon vertex can also be used [12] ; this derivation 
is, perhaps, as easy as the ghost - gluon one. 

If the Ward identities are QED-like, it is sufficient to renormalize the gluon propagator; 
no vertex calculations are required. One such case is the Coulomb gauge [2j; it has a 
slight disadvantage of being not Lorentz invariant (this makes gluon loop calculations 
more difficult). Another example is the background field formalism [13] . It is Lorentz 
invariant, and provides, probably, the shortest derivation of p in QCD. 

However, all these derivations don't explain asymptotic freedom (antiscreening) in sim- 
ple physical terms. Probably, the simplest explanation of this kind is presented in [13] (it 
is also discussed in [8]). We shall discuss in in Sect. |3j 




(1.3) 
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2 Coulomb gauge (Khriplovich 1969) 



2.1 Feynman rules 

This gauge has several advantages. It is ghost-free, and Ward identities in it are simple. 
It also has one significant disadvantage: it is not Lorentz-invariant. Currently, it is widely 
used in non-relativistic QCD [15]. In this effective theory, there is no Lorentz invariance 
from the very beginning, and the disadvantage does not matter. On the other hand, 
separating Coulomb and transverse (chromomagnetic) gluons helps to estimate various 
contributions when considering non-relativistic bound states. 
The gluon propagators in this gauge are 

^ _ = -iS ab D {q) , .www™. = -iS ab D l j (q) . (2.1) 

The Coulomb gluon propagator 

A>(?) = "4 ( 2 - 2 ) 

does not depend on q . This means that the Coulomb gluon propagates instantaneously. 
The transverse gluon propagator 

^ - ?T7o ( s " ~ f) (2 ' 3) 

describes a normal propagating massless particle (the usual denominator q 2 + iO). In ac- 
cordance with the gauge-fixing condition V- A a = 0, this 3-dimensional tensor is transverse 
to q. 

The three-gluon vertices look as usual. There is no vertex with three Coulomb legs (if 
we contract the usual three-gluon vertex with t> M = (1,0) in all three indices, it vanishes). 

We are going to calculate the potential between an infinitely heavy quark and an in- 
finitely heavy antiquark in the colour-singlet state. 

The infinitely-heavy quark propagator 



(2.4) 



Po + i0 

does not depend on p (if the mass M were large but finite, it would be 



the kinetic energy disappears at M —>■ oo). This means that the coordinate-space prop- 
agator is ~ S(f)6(t): the infinitely heavy quark does not move in space, and propagates 
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forward in time. The energy of an on-shell infinitely heavy quark is 0, and does not depend 
on p (it would be p 2 / (2M) if M were large but finite). 

An infinitely heavy quark interacts with Coulomb gluons but not with transverse ones: 



ig t a 



IWWV 



0. 



(2.5) 



The infinitely-heavy antiquark - Coulomb gluon vertex differs from (12. 5p by a minus sign. 



2.2 Quark— antiquark potential 

We calculate the scattering amplitude of an on-shell infinitely heavy quark and an on-shell 
infinitely heavy antiquark in the colour-singlet state with the momentum transfer q in 
quantum field theory (Fig. [1]), and equate it to the same amplitude in quantum mechanics. 
In the Born approximation, it is 

(2.6) 



Two-particle-reducible diagrams in quantum field theory (which can be cut into two dis- 
connected parts by cutting a quark line and an antiquark one) correspond to higher Born 
approximations in quantum mechanics, and we don't have to consider them. 

(0,0 (0,-q) 



(0,0) (0,0) 



Figure 1: On-shell scattering amplitude 
To the lowest order in g 2 , the scattering amplitude is 



-iC F glD{^q) = iC F 



9l 



(2.7) 



Therefore, 



U ? =C F g 2 D(0,q) = -C F 



9l 



(2.8) 



and we obtain the Coulomb attraction potential 



U(r) 



-Op — . 



(2.9) 



Some people prefer to discuss the quark-ant iquark potential in terms of the vacuum 
average of a Wilson loop (Fig. [2]) with T ^> r. Of course, this is exactly the same thing, 
because the infinitely-heavy quark propagator is a straight Wilson line along the 4- velocity 



v. 



Figure 2: Wilson loop 



The energy of a colour-singlet quark-ant iquark pair separated by r*is U(f). Therefore, 
neglecting boundary effects near time and T, the Wilson loop is 



-iU(r)T 



(2.10) 



or, to the first order in a s 

This order-a s term is 
T 



l-iU{r)T. 



7 ' ; ~i C F g 2 T Id (t, f)dt=-iC F gl T I - ^ ^> ' n , a , ' •/- 



(2vr) d 



(2.11) 



r 

(integration in r gives T), and hence we obtain ( 12. 9j) . 



2.3 Corrections to the scattering amplitude 

Now we want to calculate the first correction to the scattering amplitude in Fig. [TJ 

First of all, there are external-leg renormalization factors. They are given by the 
derivative of the infinitely-heavy quark self-energy at its mass shell. The infinitely heavy 
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quark only interacts with Coulomb gluons, and its self-energy is 



d^k 
k 2 



(2.12) 



(because the infinitely-heavy quark propagator does not depend on k). In the coordinate 
space, the infinitely heavy quark propagates along time, and the Coulomb gluon — along 
space. Therefore, the two vertices are at the same space-time point, and the self-energy is 



D(t = 0,f = 0) 



d d ~ l k , 



iq-r 



k- 



U(r = 0) =>• 



(2.13) 



r=0 



This is the classical self-energy of a point charge, and it is linearly divergent. In dimensional 
regularization, it is C0. 

Transverse gluons don't interact with the infinitely heavy quark, and hence there is 
only one vertex correction. It vanishes for the same reason: 



/ 

i »• 

\ 

\ 



0. 



(2.14) 



Therefore, we only need to consider vacuum polarization corrections (Fig. [3]): 

U q -=C F g 2 D(0,q), (2.15) 
where the Coulomb-gluon self-energy is 

-"#-- =-^ 2 n(g), (2.16) 

and its propagator is 



D(q) 



1 1 



<f 2 1 — 11(g) q 2 



■=(i + n( ff )) 



(2.17) 



(up to the first correction). 

^^This linear ultraviolet divergence leads to an ultraviolet renormalon singularity at Borel parameter 
u = 1/2, and hence to an ambiguity of the on-shell heavy-quark mass proportional to Aqcd [S] (see 
also ITJ], Chapter 8). 
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Figure 3: Vacuum polarization corrections 



2.4 Quark and transverse-gluon loops 

There are several contributions to the Coulomb-gluon self-energy at one loop. We begin 
with the quark- loop contribution U q . It is Lorentz- invariant, and is given by the spectral 
representation 

oo 

p q (s) ds 



q 2 — s + iO 



(2.18) 



with a positive spectral density 

p q (s)^0. (2.19) 

The propagator (and U$) is a superposition of the massless propagator and massive 
ones having various masses, with positive weights: 



Uq — —Cf^ 
q l 



+ 



-C F gl 



1 



Pqjs) ds 

q 2 + s 

p q (s) ds\ l_ 
q 2 



Pg{s) ds 
q 2 + s 



+ 



(2.20) 



Therefore, the potential U(r) is a superposition of the Coulomb potential and Yukawa ones 
having various radii, with positive weights: 



U{r) 



-a 



Airr 



1 - 



p g (s) ds 



+ / /W 



(2.21) 



The farther we are from the source, the more Yukawa potentials switch off, and the weaker 
is the interaction. We have screening. 

In QED, this is the only effect. It seems that screening follows from very general 
principles: causality (it allows one to express (I2.18P II (q 2 ) via its imaginary part p(s)) and 
unitarity (it says that this imaginary part (12. 19)) is a sum of modulus squared of transition 
amplitudes to intermediate states). The only chance not to get screening is to find some 
contribution which is not given by the spectral representation. 

Let's complete our calculation. Using Cutkosky rule (Fig. Hj) in the q rest frame, we 
see that two integrations are eliminated by two (^-functions, and 



p q (s) = T F n f 



9o s ~ 



(4yr) d / 2 V3 



+ 0(e] 



(2.22) 
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Figure 4: Quark loop discontinuity at the cut 



The ultraviolet divergence of IL is 



p q {s) ds 



s + q- 



uv 



4 ^ 9o 



s ds 



4^ « s 



Keeping only this divergent part, we obtain 



Ur, 



-c, 



q z 



1 i 4 ^ a s 
1 + 3^4^ 



+ 



[2.23) 



where dots are contributions of other diagrams. 

Qualitatively, the transverse-gluon loop is just like the quark one. Its contribution in 
the Coulomb gauge U. t (q) is not Lorentz-invariant, and the spectral representation is a 
little more complicated: 

n,(^W#^. (2.24) 
J q 2 — s + i\) 

When calculating the discontinuity by the Cutkosky rule (Fig. EJ), we cannot use the q 
rest frame, and there is one extra integration not eliminated by ^-functions. The general 
expression for pt(s, q 2 ) is rather complicated (this is the only point in which the Coulomb- 
gauge derivation is more complicated than the usual covariant one). 




Figure 5: Transverse-gluon loop discontinuity at the cut 



Fortunately, we don't need this general expression. In order to obtain the ultraviolet 
divergence of IT t (g), its limiting form at s 3> q 2 is sufficient. We can simplify the integrand 
of our single integral, and obtain 



[2.25) 
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The ultraviolet-divergent contribution of this diagram to U$ is 



u* = -c, 



9l 



1 4- l r as 

1 + 3 CA ^Te 



+ ■ ■ 



(2.26) 



2.5 Coulomb gluon 

There is one more contribution: the loop with a Coulomb gluon and a transverse one. We 
can understand its sign qualitatively. 

Let's consider infinitely heavy quark and antiquark at a distance r (in the colour- singlet 
state) and the transverse-gluon field (Fig. EJ). If we neglect the interaction, the ground-state 
energy is just 

E = U{r), (2.27) 
because the vacuum energy of the transverse-gluon field is 0. 



+ 



transverse-gluon 



vacuum 



Figure 6: Q, Q, and the transverse-gluon field 

Now let's take into account their interaction in the second order of perturbation theory. 
Transverse gluons don't interact with infinitely heavy quarks, they can only couple to 
Coulomb gluons exchanged between the quark and the antiquark (Fig. [7]). The ground- 
state energy decreases in the second order of perturbation theory. This means that the 
Coulomb attraction becomes stronger — antiscreening [T7] ! 



Figure 7: Interaction in the second order of perturbation theory 

This loop (Fig. [8]) depends on q but not on q , because we can always route the external 
energy q via the Coulomb propagator, and it does not depend on energy. Therefore, this 
loop has no cut in the go complex plane, and is not given by the spectral representation. 
Speaking more formally, we can say that the spectral density is 0, and the whole result is 
given by the subtraction term, which does not depend on q (but depends on q) . 
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Figure 8: Loop with Coulomb and transverse gluons 

This diagram can be easily calculated. Only the transverse-gluon propagator depends 
on k : 

d d k f(k,q) 



n c (g 2 ) 

The integral in /c can be taken first: 

dk 



(2ir) d k 2 + iO 



2 **o 2 -* 2 + <0 2 ^ 2 ) 



1/2 • 



We are left with a (d — 1) -dimensional integral similar to the usual d- dimensional massless 
loop. It can be reduced to T-functions via Feynman parametrization. The result is 

n^ 2 ) = C. 4 «2gC(i + 0(l)). (2.28) 

At last, let's assemble all our findings. Keeping only ultraviolet-divergent terms in the 
corrections, we have the momentum-space potential 



q 2 [ (4n) d / 2 
When expressed via the renormalized a s (fi) 



(2.29) 



this potential must be finite. Therefore, 

fa = (a - C A - ^T F n f . (2.31) 

The antiscreening term 4 comes from the Coulomb-gluon loop; it overweights the screening 
term —1/3 from the transverse-gluon loop. 

2.6 Ward identity 

Until now, we used the potential between infinitely heavy quark and antiquark. Or, if 
we cut our picture in two halves, the infinitely-heavy quark - Coulomb gluon vertex. 
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This vertex is convenient, because both the external- leg renormalization and the vertex 
corrections vanish (Sect. 12.31) . However, we can use some other vertex, for example, the 
finite-mass quark - Coulomb gluon vertex, equally easily. The tool needed to this end is the 
Ward identity. The Coulomb gauge is ghost-free. Therefore, it is natural to expect that the 
Ward identities in this gauge are simple, like in QED, in contrast to Ward-Slavnov-Taylor 
identities in covariant gauges, which are complicated by extra ghost terms. 

We shall proceed exactly as in QED. There, an external photon leg insertion with its 
polarization 4-vector parallel to its 4-momentum gives a difference of two propagators; 
most terms cancel pairwise. Now we have a Coulomb gluon, which is polarized along time. 
Therefore, let's set its incoming momentum to q — ujv, where v = (1, 0) is the 4-velocity 
of the reference frame. We denote such an external Coulomb gluon by a leg with a black 
triangle. A dot near a propagator means that its momentum is shifted by q. It is easy to 
check the identities 



CO i 



T 

U 



T 



T 

ivwvaI- — * • 



#0 



90 



T 

•- IwMM 




•UWWWI 







(2.32) 
(2.33) 
(2.34) 



where the right-hand sides are written as (colour structure) <g) (Lorentz structure) , and the 
curved arrow in (12.331) shows the order of indices in the colour structure if abc . The last 
equality (I2.34p is obvious: let's consider one of these diagrams with the transverse-gluon 
line removed; this object is a vector (has one index), and depends on two vectors, v and 
the transverse-gluon momentum p; and the transverse-gluon propagator is transverse to 
both v and p. 

In covariant gauges, the right-hand side of (12.331) contains extra terms, where one of 
the gluon lines becomes longitudinally polarized. These terms produce ghost propagators, 
thus transforming simple Ward identities to more complicated Ward-Slavnov-Taylor ones. 

Now we are ready to apply these identities to the one-loop vertex. Let's consider the 
QED-like diagram with a transverse gluon: 



UJ 



(2.35) 
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With a Coulomb gluon, 



T 

i 



S N 



— k- - ^-^ k- 



— k- - — k- 



(2.36) 



because we can route the extra momentum q along the Coulomb-gluon line instead of the 
quark one, and the Coulomb propagator does not depend on an additional momentum 
q = uv. In QED, the only contribution is (12.351) (and there is no colour factor in its 
right-hand side): the difference of the fermion self-energies with the momenta p + q and p. 
In QCD, we also have 



UJ 



go 



(2.37) 



Here we have used the definition of if abc via the commutator in order to re-write the colour 
structure as a difference: one term is the same as in (I2.35j) (and they cancel); the second 
one is the colour structure Cp of the quark self-energy times that of the elementary vertex 
it a . Two remaining contributions vanish, due to H2.34D : 



T 

i 



T 

i 



r > 







(2.38) 



Qualitatively, we can say that the non-abelian charge flows along both lines in the quark 
self-energy diagram; the longitudinal gluon insertion into one of these lines "measures" the 
charge flowing along this line; in total, the whole quark charge flows, thus giving the full 
colour structure of the self-energy Cp. 

The Ward identity provides a relation between the quark - Coulomb gluon vertex 



^ q -ig t a T(p,q), T(p,q)= l0 + A(p,q) (2.39) 

p 



at q = uv and the quark self-energy 



= -^(P) (2-40) 
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which is related to the propagator: 

S(p) 



$-m - E(p) ' 
The Ward identity can be written as 

wA(p, uw) = E(p) - E(p + uu) (2.41) 

or 

wr(p, uw) = S-^p + uu) - ^(p) . (2.42) 

We have proved it at one loop, but it stays correct at higher loops, too. 
Now let's recall how the coupling constant renormalization 

So = Ztfg 

is derived. When expressed via the renormalized coupling g, the vertex and the propagator 
should be 

r — Zy^- V j s = z^pS r , 

where the renormalized vertex and propagator are finite at e —>■ 0. The scattering amplitude 
is obtained by multiplying the proper vertex by the external-leg renormalization factors: 

gvTZ^Z 1 / 2 = gT r Z]J 2 ZyZ^Z 1 / 2 , 

1/2 1/2 

and it must be finite. Therefore, Zd ZyZ^Z a must be finite at e — > 0. But the only 
minimal renormalization constant finite at e — > is 1: 

Z a = {Z T Z, lP Y 2 Z A 1 . (2.43) 

I. e., in order to find the coupling-constant renormalization, one needs the vertex renor- 
malization factor and all external-leg renormalization factors for this vertex. 

The Ward identity makes things simpler. From f)2.42p . Z^Z^ must be finite at e — > 0, 
and hence 

ZtZ^ = 1 . (2.44) 

Therefore, 

Z a = Z~ A l . (2.45) 

The coupling constant renormalization is determined by the Coulomb-gluon propagator 
renormalization only. And this is exactly what we studied in Sect. I2.4H2.51 
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3 Chromomagnetic properties of vacuum 
(Nielsen 1981) 

3.1 Dielectric and magnetic properties of vacuum 

We shall start from QED, and then consider non-abelian theory. 

We shall use a popular approach to renormalization due to Wilson. Suppose we have 
the path integral in momentum space. We integrate out fields with momenta p > A: 

f Ud<j> p e lS = I n^ 5A , (3.1) 

J P ^ p<A 

where the effective action is defined by 

JSa = [H d<P P e tS . (3.2) 



e 



p>A 



When we are interested in processes with small momenta Pi <C A, S\ can be expressed via 
a local Lagrangian. At the leading order in 1/A, it is the standard dimension-4 Lagrangian; 
it contains renormalized fields (at the scale A) and the coupling g(A). 

Now we want to integrate out also fields with momenta between A' and A (where 
A' « A): 



js 



A' 



J] d<P p e %s \ (3.3) 

A'<p<A 

and to obtain g(A'). If we are interested in processes with characteristic momentum p, we 
should use A' not too far separated from p, in order to avoid large logarithms. 

For example, if we consider interaction between a quark and an antiquark at a distance 
r, we should use A' ~ 1/r. Then the Coulomb potential will be simply — e 2 (A')/r. But if 
we start from the theory with a high cut-off A, then vacuum modes with momenta between 
A' and A will act as a dielectric medium, and the potential will be — e 2 (A)/(er): 



e 2(A') = *^, (3.4) 



where only modes with momenta between A' and A contribute to e. If e > 1, we have 
screening; if e < 1 — antiscreening (asymptotic freedom). 

In the case of an ordinary matter, its dielectric and magnetic properties are independent. 
But vacuum should be Lorentz- invariant. Signals should propagate with velocity 1: 

efx = l. (3.5) 

Therefore, diamagnetic vacuum (// < 1) means screening, and paramagnetic one (// > 1) 
— asymptotic freedom. 
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When we switch magnetic field B on, the vacuum energy changes by 



We shall show that 



This means 



and therefore 



AE vac = 



2 A 2 



\i = e 1 = 1 + A 
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e 2 (A'' 



(4tt) 2 ° A' 2 
A 2 



1 + /5o (47F log A' 2 



e*(A) 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



i.e., /?o is indeed the 1-loop /5-function coefficient. 

The vacuum energy of a charged scalar field (describing particles and antiparticles) is 



EL 



(3.10) 



For a charged fermion field, it is the energy of the Dirac sea 

ffvac = - y~]wj. (3.11) 

i 

Therefore, in general we can use ( 13. 10D multiplied by (— l) 2s . 
3.2 Pauli paramagnetism 

How the vacuum energy changes when we switch the magnetic field B on? There are two 
effects: spin and orbital, and they can be considered separately. 

First we discuss interaction of the spin magnetic moment with the magnetic field. 
Without the field, a massless particle has the energy uo = k. When the magnetic field is 
switched on, the energy becomes 

u = y/k?- g s s z eB , (3.12) 

where g s is the gyromagnetic ratio for our spin-s particle (we shall discuss this in Sect. 13.61 
in more detail). 

Suppose the magnetic field is along the z axis. Massless particles only have 2 spin 
projections s z = ±s. The vacuum energy change at the order B 2 is 

Vd 3 k 



AE- 



Pauli 



\2s 



-1) 2S V 



(2tt) 3 
{QsSeBf 



\/k 2 + g s seB + \/k 2 - g s seB - 2k 



-2(-lY s (g s s) 



(47T) 



d 3 k 1 

J2^¥ ~ ~ ( 
i A2 B \r 



T (3 - 13) 
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where only modes with momenta between A' and A are included. Let's stress once more 
that what we are calculating is the vacuum energy: there are no particles, only empty 
modes. 



3.3 Landau levels 

Now let's discuss the effect of magnetic field B on the orbital motion. In order not to have 
complications related to spin, we consider a massless charged scalar field (p. Its Lagrangian 
is 

L = {D^) + D^. (3.14) 

For magnetic field B along the z axis, we can choose the vector potential as A y = Bx, 
A x = A z = 0. Then the equation of motion is 



V 2 - e 2 B 2 x 2 - 2ieBx^- + E 2 



Its solutions have the form 



(f = e 



dy 



if = 0. 



(3.15) 



(3.16) 



The equation for <p(x) has the same form as the Schrodinger equation for harmonic oscil- 
lator: 



i d 2 - 2 



2dx 2 + 2 n 



(3.17) 



The oscillator energies are 

E n = uj{n + \). (3.18) 

Comparing A3. 1 5[) with f)3.17p . we see that the energies of our massless particle in magnetic 
field B are 

E 2 = k 2 z + 2eB (n + |) , (3.19) 
and the corresponding wave functions are 



if = e 



i(k y y+k z z) 



(3.20) 



In other words, E 2 consists of discrete Landau levels of transverse motion plus free motion 
along the magnetic field. 

Each Landau level has a high degree of degeneracy. In order to find it, let's put our 
particle into a large box V = L x X Ly X L z . Then the allowed longitudinal momenta are 



therefore, the number of allowed modes in the interval dk 7 is 



drt, 



Lzdh 
2tt 
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Similarly, the allowed values of k y are 



ky Tly . 

L,y 



As we see from (I3.20p . k y is related to the x coordinate of the center of the Larmor orbit. 
It must be inside our box: 

s 6 ^ 



Therefore, 



n y G 







eB L x L r 
2tt 



Energy does not depend on k y . Hence the degeneracy of each Landau level is 



eB L x L y 
2^ 



(3.21; 



It is equal to the magnetic flux through our box (B L x L y ) measured in flux quanta 27r/e. 

The spectrum of E\ = E 2 — k 2 z at B = is continuous. The number of states in the 
interval dE\ is 

LxLy -dE 2 . 



4tt 



When the magnetic field B is switched on, each interval AEj_ = 2eB is contracted into 
a single Landau level (Fig. [9]). The number of states in each interval is the Landau level 
degeneracy (|3.21|) . 



AEl/{2eB) 




Figure 9: Continuous spectrum at B = and Landau levels 
The vacuum energy of our massless charged scalar field in the magnetic field B is 



(3.22) 



n=0 
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where 



/(*) 



eBV 



+00 



^k 2 z +2eBxdk z 



(3.23) 



3.4 Euler summation formula 

Integrals are more convenient than sums. If we have a smooth function f(x) (i.e., its 
characteristic length is L 3> 1), then, obviously, 



N N+1 

/ /( 

n=0 J n 



x) dx . 



But how to find a correction to this formula? 

Let's re-write this integral as a sum of integrals over unit intervals: 



(3.24) 



N+1 N 1/2 

I = J f(x) dx = J f (n+ ~ + x) dx . 



n_0 -l/2 

The smooth function f(x) can be expanded in Taylor series in each interval: 

1/2 



N 



n=0 



f(n + \) + -f(n + \)x> + . 



dx 



1/2 



(terms with odd powers of x don't contribute). Calculating the integrals, we get 



N 1 N 



n=0 



n=0 



The second sum in the right-hand side is a small correction (because /" ~ //L 2 ); therefore, 
we can replace it by an integral (see (I3.24I) ): 



N N+1 N 

^ = E/(^+D+^ / f"{x)dx + ... = Y^f{n+\) + ^f{x) 

n=0 n n=0 



N+l 



+ 



Finally, we arrive at the Euler summation formula 



N N+1 

J2f( n + D= I /(*) dx - -f(x) 

n=0 i 



N+1 



+ 



(3.25) 



The correction here is of order 1/L 2 . It is easy to find a few more corrections, if desired. 
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3.5 Landau diamagnetism 

The vacuum energy in magnetic field B ( 13.221) can be re-written using the Euler for- 
mula (13.251) as 



£vac = J f(x) - Lf( x ) 




(3.26) 



The integral here is the vacuum energy at B = (when the spectrum of E\ is continuous, 
Fig. [H]). The shift of the vacuum energy due to the magnetic field is 

1 w,^ e 2 B 2 V fdk 1 e 2 , A 2 B 2 , n _ 

where only modes with momenta between A' and A are included^]. 
3.6 The QED result 

The full AE vac in QED is the sum of the spin contribution (13. 131) and the orbital one (13.271) : 

A£ vac = -/3 07 ^log^-!V, (3.28) 

where 

(3 = 2J2(-I) 2s [(9ss) 2 -^-} (3.29) 

s 

is the sum over all charged fields, and n s is the number of polarization states: no = 1, 
n s ^ = 2. 

Let's demonstrate that a Dirac particle (e.g., electron) has gyromagnetic ratio gi/2 = 2. 
For a massless particle, the Dirac equation is 

Pi) = , D„ = - icA^ . (3.30) 

Let's multiply it by Jp: 

l/) 2 ij = 0, Ip 2 = d 2 - e 2 A 2 - 2ieA"d fl - ie-f-fd^ 
(in the last term, acts only on A„). If we suppose that d ■ A = 0, 

We choose A^ = (0, 0, Bx 1 , 0), then F 12 = -F 21 = -B and 

]jp = D 2 + ieBj 1 -/ 2 = D 2 + 2eBs z . 



2 Many subtleties have been swept under the carpet in this derivation sketch. A somewhat more accurate 
derivation [14] yields the same result. 
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The equation of motion in magnetic field B (directed along z) becomes 



d 

V 2 - e 2 B 2 x 2 - 2ieBx— + 2eBs z + E 2 

dy 



if) = 



Its energy spectrum is 

E 2 = k 2 z + 2eB (n + ±) - 2eBs z 
(see Sect. !3.3jl R Comparing this with (13.121) . we see that 



01/2 



2. 



(3.31) 
(3.32) 
(3.33) 



So, the electron contribution to /3q (I3.29|) is —4/3; if a charged scalar particle exists, 
it contributes —1/3. Note that for s — 1/2 the spin effect overweights the orbital one. 



However, due to the factor 



\2* 



the spin effect leads to diamagnetism, and the orbital 



one to paramagnetism. This is because here we are interested in the Dirac see. In physics 
of metals, we are interested in positive-energy electrons (below the Fermi surface), and the 
spin effect gives Pauli paramagnetism, while the orbital one — Landau diamagnetism. 



3.7 The QCD result 

Now it is easy to obtain (3 in QCD. We have chromomagnetic field instead of magnetic. 
Let's choose its colour orientation along an axis a such that t a ° is diagonal (for the SU(3) 
colour group with the standard choice of the generators t a , t s is diagonal, and we choose 
ao = 8). 

The quark contribution follows from the electron one in QED. The contribution to /3q 
is proportional to the charge squared. The sum of squares of colour "charges" of a quark 
is Tr t a °t a ° (no summation). Recalling 

Tr t a t b = T F 5 ab , 

we arrive at the contribution 

- (1 - |) 2T F n f (3.34) 
of 71/ quark flavours. If there were scalar quarks, each flavour would contribute 




What about gluons? First of all, we have to find their gyromagnetic ratio g±. We shall 
do this for the SU{2) colour group, because calculations are simpler in this case; the result 
will be valid for other colour groups, too. Let's consider the SU{2) Yang-Mills equation 

D U G% = (d u S ab + ge acb A cu ) G\ v = . (3.35) 

3 Strictly speaking, we should consider spin and orbital effects together, using (|3.32[) . It is easy to see 
that in the order B 2 they can be treated separately. 
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The external field is A? . We linearize in the small components A 1 ^ 2 . For A^ = A 1 ^ — iA 2 ^ 
we get 

D v G^, + igG%,A- v = 0, 

where 

G> = D^A~ - D V A~ , D„ = d lx - igAl . 
In the D^A~ = gauge, the equation of motion becomes 

D 2 A- - <ligG\ v A- v = 

(we have used [D^,D U ] = —igG 3 ^). Our external field is oriented along z in space and 
along 3 in colour: G\ 2 = — G\ x = —B. Using A~ 2 = iSzA -1 (s z = ±1), we finally obtain 

[D 2 + 2igBs z ]A- 1 = 0. (3.36) 

This equation looks exactly the same as (13.311) . and hence g\ = 2. 

Gluons with colour a\ such that t ai is diagonal don't interact with our chromomagnetic 
field (for the standard SU(3) generators, a x = 3). All other gluons can be arranged into 
pairs with positive and negative "colour charges" (particles and antiparticles). The sum 
of their squares (both signs!) is Ca- in the adjoint representation 

Tr t a t b = C A S ab , 

i.e. we have to replace 2e 2 — > C^g 2 in QED results. 
Finally, we arrive at 

0o = (4 - |) C A - (1 - |) 2T F n f . (3.37) 

Pauli paramagnetism of the gluon vacuum (g\ ■ l) 2 = 4 is stronger than its Landau dia- 
magnetism —1/3. This leads to antiscreening (asymptotic freedom). 
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